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Planetary Isophotes as a Clue to Aerosol Characteristics 
By George W. Kattawar and Andrew T. Young 

ABSTRACT 

A study was made to see how much information could be extracted from the 
Mariner 9 Mars isophotes taken at a phase angle of -60°. We found that the 
Minnaert functions and both isotropic and Rayleigh scattering could easily be ruled 
out, and that it was essential to use forward-peaked phase functions, which were 
computed from Mie theory. We could achieve isophotes similar to those observed 
assuming a semi-infinite dust cloud with a considerable variation in particle 
properties and size distribution, so long as the ratio of the multiply- to 
singly-scattered light was held within certain limits. These conditions are 
met by micron-sized, moderately absorbing mineral grains whose mean size 
should not be much larger than a micron. It was also found that a dust, cloud 
of finite optical thickness bounded from below by a Lambert ground would fit 
the isophote data. 


Planetry lsophotes as a Clue to Aerosol Characteristics 
I. Introduction 

Spacecraft television systems can, in principle, provide photometric infor- 
mation at much better spatial resolution than ground-based telescopes, v/hich are 
severely hindered by atmospheric blur (typically, 0.1 planetary radius or worse.) 
However, past television systems have failed to produce reliable photometry: 
reduced data were found not to be proportional to scene brightness (see, e.g., 
Young, 1974a.) In spite of this quantitative failure, some qualitative infor- 
mation remains. In particular, reduced data numbers do seem to have a smooth 
monotonic relation to scene brightness, so that the shapes of isophotal contours 
in the TV data are approximately correct, even though the numerical values of 
each contour level are quite uncertain. 

Such contours have been published from Mariner 9 data (Young, 1974a, Fig. 9) 
and Venus from Mariner 10 data (Devaux et al , 1975, Fig. 2). In this paper, we 
shall concentrate on interpreting the Mars isophotes. 

II. The Data 

Fig. 1 shows the adopted contours in a standard format. We assume the 
north-south symmetry required by the reciprocity principle, and show contours 
for only one hemisphere. In this and subsequent figures, the points C,L,M,P,S, 
and T are the disk center, the limb, the mirror point the photometric pole, 
the subsolar point, and the terminator, respective. The phase angle for 
the Mars data is approximately 60°. 

We may summarize the principal features of these isophotes as follows: 

1. They are all convex curves; their limb sides are nearly parallel to the limb, 
and their terminator sides all curve the same way as the terminator. 
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2. Each isophote has a rather sharp corner at the junction of its limb and 
terminator portions. 3. The terminator side of each isophote is much 
straighter than the limb side, so that they are shaped like a letter D. This 
is particularly true of the isophotes between the points C and M. 

As pointed out by Young (1974a), the position of the maximum brightness 
is very near the subsolar' point, S, which is characteristic of a Lambert photo- 
metric function. However, the close adherence of the isophotes to the limb, 
up to high latitudes, suggests a large single-scattering contribution in the 
reflected light. If one were to try to use a Minnaert function to describe the 
brightness distribution, a value of k substantially below unity would be re- 
quired to produce this limb-brightening effect. But this would displace the 
peak brightness strongly toward the limb from the subsolar point (see Fig. 2). 

A Minnaert function with k < 1 is also unphysical, because it produces an in- 
finite brightness at the limb. Therefore a more realistic, physical model is 
required. 

A realistic scattering model involves many parameters. At the very least 
we need an anisotropic phase function, characterized by the anisotropy and the 
single-scattering albedo; if a Mie-theory model is used, these are related to 
the size distribution and the complex index of refraction of the aerosol parti- 
cles. The size distribution involves at least two parameters (mode size and 
standard deviation), and so does the refractive index (real and imaginary parts.) 
It is not clear that so many parameters can be extracted from our data, so we 
must use arty other information available. 

It would be very helpful if we knew the absolute level of each isophote, 
as well as the Bond albedo of the planet at the time of the dust storm. 
Unfortunately, as Young (1974a) points out, not even the relative intensities of 
these isophotes are accurately known. Nor is there any accurate estimate of 
the planet's Bond albedo at the time of the dust storm. However, we can 
approximate the latter as* follows: Young (1974b) has analyzed UBV ground-based 
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observations of Mars and has shown that a substantial brightening occurs at 
L s ~270°, during the annual dust storm. The average brightening Is about 0.3 mag 
in U, over 0.2 mag in B, and nearly 0.2 mag in V. As the 1971 dust storm was 
unusually intense* we can expect at least 0.2 mag brightening, and perhaps con- 
siderably more, at the time Fig. 1 was obtained. 

The average Bond albedo of Mars is about 0.08 in B, 0.15 in V, and 
o 

0,25 at 6264A (Irvine et al., 1368). If we knew the effective wavelength of 

Fig. 1, we could estimate an average Bond albedo for the planet at this 

wavelength. According to Young (1974a) the effective wavelength of the B camera 

o 

for sun! ig ht is about 5590A. Mars is much redder than the Sun, so the effective 

wavelength will be greater than this. The average (B~V) color of Mars is about 

2/3 of the way from that of a K5 V star to that of a MO III star; according to 

Azusienis and Straizys (1966), the effective V wavelengths for these stars are 
o o 

5540 and 5564A, respectively , while that for the Sun (G2 V) is 5490A. According 

o o 

to these figures, the V effective wavelength for Mars would be about 5556A, or 66A 

redder than for the Sun. Since the B camera's effective wavelength for sunlight 
o o 

is 5590A, we would expect (roughly) a 66A greater effective wavelength for Mars, 

• o 

or 5656A. The interpolated Bond albedo at this wavelength is about 0.16. 

On the other hand, Leovy et al., (1972) give an effective wavelength for 

o 

(Mars + the Mariner 9 B Camera) of 5850A, at which the albedo would be about 
0.i8. There Is considerable uncertainty in the interpolation, because of the 
nonlinear run of the Martian albedo with wavelength. We adopt 0.17 ± .02 as a 
reasonable value for the Bond albedo of Mars as seen by the B camera, withou t a 
dust storm. The storm probably increased the Bond albedo by at least 0.03, so a 
minimum value of the Bond albedo, = 0,20 is appropriate for our data. Values 
as large as 0.25 may not be excluded, considering the uncertainties. 


II Method of Computation 


Referring to Fig. 3, let OE and OS be the direction to the satellite and 
Sun respectively; SEO the plane defining the Intensity equator and also the 
phase angle a, Q an arbitrary point on the planetary surface where OQZ defines 
the vertical at point Q. . It Is also convenient to introduce the planetary co- 
ordinates ^(longitude) and n(colatltude) . 

The angles of Incidence and reflection can be expressed as 

cos 6 q - p q » sin n cos(?-a) (1) 

cos o * u * sin n cos r Jt (2) 

where 0 Q and e are acute angles. 

The azimuthal angle <{> for the reflected ray is given by 

cos (cos 0 o - cos «)/sin 0 sin f) Q (3) 

and 

sin *= sin “cos nsin 0/sin 0 Q . (4) 

The x and y coordinates of the projection of point Q on the plane OPG are given 
by 

x = sin n sin 5 (5) 

and 

y = cos n , ( 6 ) 

where we have assumed a planetary radius of unity. 

To generate isophotes, we used a grid in x and y coordinates, related to 
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* n 11 0 8 J1 cos a + x sin a. (8) 

Once y, y Q , and $ have been determined, the reflected intensity I (y, y Q , <?>) can 
be found for the point (x,y). The intensities were computed by the matrix 
operator routine described by Plass et al., ( 1973 ). To insure adequate reso- 
lution in both y and <{>» twenty-one Lobatto points were used for the y coordinate, 
and nineteen equally-spaced azimuthal angles (a$ = 10°) were used for the 4> coordi 
nate in the computation of the reflected intensities. 

As an aid in interpreting the data, isolines of y, y Q and 4. are presented 
in Fig. 4 . They were computed by the following procedure. Rewriting Eq. ( 7 ) 
we get 

x 2 + y 2 ** 1 - y 2 • ( 9 ) 

This shows that the isolines of y are circles, centered at the center of the 
disk and having a radius of (1 - y 2 )^ 2 . Now eliminating y from Hqs. (7) and 
(8) we get 

(x - y Q sin a ) 2 y 2 = (1 - y Q 2 ). 00) 

COS 2 a 

This shows that the isolines of y Q are ellipses whose centers are located at 

(y Q sin a, 0 ). The ratio of the major to minor axes of these ellipses is 

(cos a)" 1 . For the phase angle of 60° shown in Fig. 4 this ratio is simply 

2:1. It should be noted that the isolines of u Q are also isophotes for a Lambert 
surface since the intensity for a Lambert surface is proportional to y . The 
isolines of 4' are more difficult to compute and the following technique was 
used. We have: 

tan o/2 = cos [(<j> - 3 )/] tan £(ct - o o )/2]/cos £(4+3 )/2j (11) 


t 


and 


B 


sin 0 O * sin 0 sin o/sin 4 . ( 12 ) 

First the angles 4 and 0 were chosen, and Eq. (12) was used to compute 0 O . 

These values were then used in Eq. (11) to compute the angle 0 . Once 0 and 
e Q have been computed then the x and y coordinates can easily be calculated 
from Eqs. (7) and (8). A similar technique was used by SeKera and Viezee (1961). 

The phase functions used in this study were generated from Mie theory with 
a program developed by Kattawar and Plass (1967), using a particle size distri- 
bution of the form 


n(r)« r a exp (-br Y ), 


03) 


which Deirmendjian (1969) called a modified gamma distribution. The mode radius is 

1/Y 

r c ■ ( fb> • < 14 > 

The spherical or Bond albedo Ag can easily be computed from the reflected 
flux F R (n 0 ), for an incident solar flux of unity, by 

" 2 V^o^V 

Now if one knows the Bond albedo for a certain phase function with single 
scattering albedo iXj, then one can use the similarity relations of Van de Hulst 
(1968) to relate it to another phase function with single scattering albedo ^ 
by 


u 


1 _ 


(ii. 


(16) 


1-9! 1 - g 2 

where g, the so-called asymmetry factor, is related to the phase function 
P(cos 0 ) by g = P(cos 0) cos 0 d(cos 0 )/ ^ P(cos 0 ) d(cos 0 ). (17) 


It should be noted that Eq. (16) is only approximate, and the error becomes worse 
as w decreases. When to is close to unity it is quite accurate. 


IV. Discussion of Models 
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The model Isophotes are presented in a uniform format. The contour In- 
terval is 0.08 of the maximum brightness of each model; the relative brightness 
levels are thus 0.08, 0.16..., 0.88, and 0.96. The 0.48 contour, which has 
half the brightness of the brightest contour, is emphasized. To show more de- 
tail near the peak, the 0.92 contour level is usually added as a dashed line. 

Thus the 12 solid and 1 dashed contours are comparable to the 14 solid contours 
in Fig. 1. Only one hemisphere is shown; the other would be its mirror image. 

Fig. 2 compares several very simple models. The top row shows isophotes 
for Minnaert functions; the middle shows isotropic-scattering models; and the 
bottom row shows Rayleigh-scattering models. As mentioned above, the Minnaert 
functions fail to reproduce the data. A further problem is that a Minnaert 
function has no definite albedo, so this information cannot be used. 

Isotropic-scattering models of infinite optical depth were proposed by 

Leovy, et al, (1972). According to them, a single-scattering albedo i^of about 

0.7 can reproduce the surface brightness measured by the B camera. However, a 

subsequent investigation of the TV photometry (Young, 1974a J casts doubt on the 

reliability of the data on which this result was based. Fig. 2f shows that the 

shape of the isophotes calculated for a semi -infinite isotropic-scattering model 

with W Q = 0.7 bears no resemblance whatever to the observed shapes, which are 

unaffected by the nonlinearity of the data. However, the Bond albedo for this 

* 

model (0.26) is acceDtable. The conclusion of Leovy et al . (1972) that isotropic 
scattering produces a quasi -Lambertian photometric function for such a low 


*The Bond albedo of the ui 0 = 0.7 model is about 50f' brighter than the albedo with- 
out a dust storm. Thus, the albedo derived for the dust by Leovy et al . (1972) is 
inconsistent with their statement that "the dusty Mars is not significantly 
brighter than the clear Mars." 


0 
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single-scattering albedo Is clearly wrong. Fig. 2d shows that quasl-Lambertlan 
behavior does occur for conservative Isotropic scattering; however, even Fig. 2e, 
with a Bond albedo roughly twice th'.t of Mars, is far from Lambertian. 

These isotropic models already show the advantages of using Isophotes, 
rather than the unreliable Mariner photometry, as a basis for studying the 
optical properties of the dust. Clearly, no Isotropic-scattering model will fit 
the data. 

Another simple (but more plausible) model is Rayleigh scattering. Although 

„3 

the optical depth of the molecular atmosphere of Mars is only about 3 x 10 , 

very small dust particles could produce nonconservative Rayleigh scattering. 
However, as the single-scattering phase function is not very different from iso- 
tropic, we do not expect the computed isophotes to be much like those observed. 
Figs. 2 g-i confirm this expectation. Conservative Rayleigh scattering, like 
conservative isotropic scattering, produces quasi-Lambertian Isophotes; but 
models with plausible Bond albedos are much too bright at the limb. 

Both of the conservative-scattering models (Figs. 2d and 2g) have their 
maximum brightness closer to the mirror-point M than for a Lambert surface, 
whose maximum is at S, the sub-solar point. In these models, most of the ob- 
served light (which is quasi-Lambertian in distribution) is multiply-scattered. 

The ratio of multiple- to single-scattered intensity for the brightest point 
(MS/SS) is 6.7 for w 0 = 1.0, whereas for us Q = 0.80 and 0.5710 it is 1.8 and 
1.4 respectively. Hence, increasing the importance of single scattering should 
move the maximum to S and produce limb brightening. Not only do the observed 
isophotes indicate a moderate brightening toward the limb, but the non-conservative 
Isotropic and Rayleigh models of low Bond albedo show strong limb brightening, 
due to a much smaller multiple- to single-scattering ratio. We need a low-albedo 
model that has relatively less single-scattering backward at 60° phase angle. 

Any forward-scattering phase function should do this. 


In order to use Mte-scatterlng models we have to choose both the com- 
plex refractive Index and the size distribution. In Figs. 5a - 5c the 
size distribution Is fixed and the refractive index Is varied. The phase 
functions for these models are shown in Fig. 6a. In Fig. 5a the Bond albedo 
Is too high and the Isophotes are too rounded, due to the large value of 
MS/SS ■ 36.2. To decrease the Bond albedo, we Increased N g (Fig. 5b) and 
thus decreased w Q and reduced the MS/SS ratio to 16.3. This reduction 
Is due entirely to the diminution of the multiple-scattering, since the 
single scattering contribution Is approximately the same for the two models 
(see Fig. 6a). Consequently, the brightest isophote moves closer to the 
limb arid has a more realistic D shape, while the Bond albedo decreases to 
the estimated value during the dust storm. In Fig. 5c we increased the real 
part of the refractive index, and adjusted the imaginary part to keep the Bond 
albedo in agreeing., with r ig. 5b. The isophotes are shifted more toward the 
limb and extend higher in the vertical direction; also, the M3/S5 ratio has 
been reduced to 8.5, which is largely due to the increase in the single scatter- 
ing contribution (see Fig. 6a). In Figs. 5e and 5f we set the real part of 
the refractive index to 1.5 and varied the size distribution. The mode 
radii for the two cases are 0.5 urn and 0.2 pm r pectively (see Fig. Ga for 
the phase function). For both of these models the isophotes are too rounded 
at the top. The smaller sizes also shift the bright point towards the terminator. 

It should also be noted that using the similarity relations of Van do 
Hulst (1968) can lead to erroneous results. For example, in Figs. 5e and 5f 
the ratio (1 - w Q )/(l - g) for the two cases is 0.25 and 0.24 respectively. 

This relation is useful for matching Bond albedo data, but is useless for iso- 
phote matching; this conclusion was already reached by Kattawar (1975). An ex- 
treme example can be seen by comparing Figs. 5e, 5f and 2h. f r the Rayleigh 
case (1 - w Q )/(l - g) = 0.2 and gives the same Bond albedo, but the isophotes 
bear no resemblance to those in Figs. 5e and 5f. 


, which has a mode radius of 


fi 1 JV 

We next used a size distribution of re 
0.5 mi, and Increased the real part of the refractive Index to 2.0. This 
has the effect of increasing the single scattering at 120° scattering angle 
(compare Figs. 6a and 6b). This gives the same Bond albedo as that of Fig. 5e 
with the same distribution but 1.5. The effect of this (see Fig. 7c) is 
to shift the maximum toward the limb and produce less rounding of the isophotes 
at higher latitudes. To bring the Bond albedo down to a more realistic level 
we increased N 2 (see Fig, 7d). The MS/SS ratio is reduced and the isophotes 
are shifted even more toward the limb and are actually starting to become con- 
cave. 

In all of the models we have considered to this point, we have assumed the 
optical thickness to be infinite. Although the observed isophotes in the sub- 
solar area are quite uncertain, both because of the low brightness gradients and 
the large residual -image effects, we suspect that effects of finite optical 
thicKuess may appear there, as the total air mass is least at the mirror point. 
We therefore explored the effect of a finite optical thickness on the shape of 
the isophotes. We used the model in Fig. 7d and first sett := 4 with a 
perfectly absorbing ground (i.e., the ground albedo = 0.0.) The results are 

shown in Fig. Qa. The Bond albedo is in an acceptable range, but the isophotes 
are pushed much too close to the limb, and have become concave over a region 
from the mirror-point to the sub-solar point, ,'e next set Aq = 0.5, assuming 
a Lambert ground, and the Bond albedo increased only by 6*'. However, the change 
in the isophotes was remarkable (see Fig. 8b). In fact, this case gave reason- 
ably good agreement with the actual isophotes. We considered two additional 
cases, namely Aq = 1.0 and a perfectly antispecular ground (see Figs. fle and 
ud). In each of these cases the isophotes are too rounded. 
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Conclusion 

The isophotes mainly tell us the ratio of multiple to single scattered 
light at our phase angle. We can achieve isophotes similar to those observed 
with a considerable variation in particle properties and size distribution, 
so long as this ratio Is achieved, and the Bond albedo is not unreasonable. 
These properties are satisfied by micron-sized, moderately-absorbing mineral 
grains. We b?lieve that the mean particle size should not be much, larger than 
a micron, for this would make too much forward scattering and put the maximum 
brightness on the limb side of the sub-solar point. 
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Figure Captions 


Figure 1. Mariner 9 Isophotes of Mars during dust storm, at a phase angle 
-60°. 


Figure 2. 


Figure 3. 
Figure 4. 

Figure 5. 


Figure 6. 


Figure 7. 


Figure 8. 


a-c: Computed- Isophotes for Mlnnaert functions, where the intensity 

k k-1 

Is assumed proportional to ^ y for k - 0.9, 0.8, and 0.7, 
respectively, d-f: Isophotes for a semi - Infinite isotropic scattering 
atmosphere with single scattering albedos = 1.0, 0.9, and 0.7 
respectively: g-i: Isophotes for a semi-infinite Rayleigh scattering 

atmosphere with = 1.0, 0.8, and 0.57 respectively. 

Planetary coordinates used for computation. 

Isolines of y, y o , and <j> for a phase angle of 60°. The isolines 
of y Q are isophotes for a Lambert surface. 

Isophotes for various Mie scattering models assuming a semi- 
infinite cloud, a-c: n(r) re r®e~^ r and N-j - 1.5, Ng - 0.001; 

N-j « 1.5, N 2 = G.004; and Nj «* 1.65, N 2 = 0.005 respectively, 
d: Mars data, e: n(r)<* r 6 e"' 12r , N-j = 1.5, N 2 = 0.005. 
f: n (r ) “ r 6 e“ 30r # .N 1 « 1.5, Ng « 0.013. 

a-b: Phase functions computed from Mie theory for use in isupiiote 

computations. 

a-b: Same as Figs.. 5e and 5d respectively, c-d: n(r)«r G e”^ 2t and 

N-| = 2.0; N 2 = 0.0^55, N 2 = 0.01 respectively. 

Isophotes showing the effects of a finite optical depth,!: = 4, over 
both a Lambert surface (a-c) and an antispecular surface, (d ) . 
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